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ABSTRACT 

Carbon-oxygen white dwarfs contain ^^Ne formed from a-captures onto '^N during core He burning in the 
progenitor star. In a white dwarf (type la) supernova, the ^^Ne abundance determines, in part, the neutron- 
to-proton ratio and hence the abundance of radioactive ^^Ni that powers the lightcurve. The ^^Ne abundance 
also changes the burning rate and hence the laminar flame speed. We tabulate the flame speedup for different 
initial '^C and ^^Ne abundances and for a range of densities. This increase in the laminar flame speed — about 
30% for a ^^Ne mass fraction of 6% — affects the deflagration just after ignition near the center of the white 
dwarf, where the laminar speed of the flame dominates over the buoyant rise, and in regions of lower density 
^ 10^ gcm^^ where a transition to distributed burning is conjectured to occur. The increase in flame speed will 
decrease the density of any transition to distributed burning. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — supernovae; general — white dwarfs — 
galaxies: evolution 



1. INTRODUCTION 

In the past decade, type la supernovae (hereafter, SNe la) 
have become the premier standard candle for measuring the 
geometry of the universe. Although many details of the 
explosion are not well understood, there is a general be- 
lief that the explosion is the thermonuclear incineration of 
a C/O white dwarf that has increased in mass through ac- 
cre tion to j ust below the Chandrasekhar limit (for a review, 
see iHillebrandt & Niemeyer 2000) . Current obse rvations are 
samp ling the SNe la population out to z w 1.6 dRiess et alj 
120041) . and future missions will push this limit even farther to 
z < 2. The larger sample of SNe la carries with it the prospect 
for discovering the progenitors of these events and their evo- 
lution towards igniti on. Numerical models (for a sampling 
of recent work, see iGamezo et all l2004t iPlewa et alJ 120041: 
iRopke et alJ l2006 l) are steadily becoming more refined and 
can begin to probe the connection between the properties of 
the progenitor white dwarf — its birth mass, composition, and 
binary companion — and the outcome of the explosion. 

The composition of the progenitor white dwarf should 
play a role in setting the peak brightness and the composi- 
tion of the ejecta. The C:0 ratio is se t by the mass of the 
progenitor main-seqence star, although Ropke & Hillebrandt 
12004.) find that the C:0 ratio is of secondary importance 
in setting the explosion energetics. After '^C and '^O, the 
next most abundant nuclide is ^^Ne, which is synthesized 
via ^'*N(a, 7) '^F(j3+)'*^0(a, 7)^^Ne during core He burning. 
The abundance of ^^Ne is therefore proportional to the ini- 
tial CNO abundance of the progenitor main sequence star 
iTimmes et alJ (l2003l) showed that the mass of ^^Ni synthe- 
sized depends linearly on the abundance of ^^Ne at densities 
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where electron capture rates are much slower than the explo- 
sion ti mescale, ~ 1 s. Simulations with embedded tracer par- 
ticles dTravadio etal J 120041: iBrown etal J 120051: iRopke et alj 
|2006|) have confirmed this dependance. 

These simulations studied the effect of adding ^^Ne by 
post-processing the (p,r) traces, and as a result did not ac- 
count for variations in either the progenitor structure or the 
sub-grid flame model caused by changes in the ^^Ne abun- 
dance. One-dimensional studies that did atte mpt to incor- 
porate diff erent progenitors self-consistently (IHoflich et alJ 
[l998 : Domfnguez et al.ll2001h found a much smaller depen- 
dence of the ^^Ni yield on metallicity^. 

The possibility that type la supernovae might evolve with 
the abundance of a-elements in the host population, com- 
bined with questions about whether this introduces system- 
atic variations in the Phillips relation, motivates further in- 
vestigation of how the progenitor composition influences the 
explosion. As a first step, we investigate in this Letter how 
the abundance of ^^Ne affects the laminar flame speed S'lam 
and width 5iam of a '^C-^^O-^^Ne mixture. Our principal 
conclusion is that S\am increases roughly linearly with the 
^^Ne mass fraction X22. At X22 = 0.06, the speedup varies, 
for carbon mass fraction X\2 — 0.5, from « 30% at densities 
> 5.0 X 10*^ gcm"^ to « 60% at lower densities. 

These calculations are relevant for two regimes: 1) the 
initial burn near the center of the white dwarf where the 
gravitational acceleration is small and the laminar flame 
speed dominates the evolution of a bubble of ignited ma- 
terial (see, e.g., Zingale & Dursi 2006), and 2) the burning 
at densities 10^ gcm^^ where the Gibson length scale 
becomes Iq ~ 5iam- The Gibson scale Iq is defined by 
v(^g) — Su m where v(£) is the eddy velo city for a length- 
scale I (see IHillebrandt & Niemeverll2()0 0, for a succinct re- 
view). The region where 5iam = is conjectured to be a 

^ It is unclear whether these studies allowed [O /Fe] to vary as a function 
of [Fe/H]. 



2 



Chamulak et al. 



TABLE 1 
430-NucLiDE Reaction Network 
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possible location for a defl agration-to-detonation transition 
jNiemever & WooslevUlQQTl) . Our calculation does not apply 
in the flamelet regime, where the buoyancy of the hot ashes 
generates turbulence via the Rayleigh-Taylor instability. In 
this regime, the effective fr ont speed becomes independent 
of the laminar flame spee d (iKhokhlovll 1 995t iReinecke et alJ 
Il999t IZhang et aP l2006l) . and the composition affects the 
front speed only through the Atwood number. At = (pfuei — 
Pash)/(Pfuei +Pash), where Pfuei(ash) is the density in the un- 
burned (burned) material. 

In § |2] we describe our computational method and 
benchmark our calcula tions against earUer results of 
iTimmes & WooslevI (11992 ). Section |3]presents the computed 
flame speeds as functions of Pfueb and X22- We provide 
a fitted expression for 5iam as a function of these parameters. 
We also give a physical explanation for the speedup before 
concluding, in § |4] with a discussion of how the transition to 
distributed burning would occur at a lower density if the ^^Ne 
abundance were increased. 

2. THE LAMINAR FLAME 

To solve for the conductive flame speed 5iani, we used the 
assumption of isobaric conditions to cast the equation for the 
energy as two coupled equations for the te mperature and flux 
(ITimmes & Wooslevl[T992l:lBildsten|[T995l) . 
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Here F is the heat flux and Cp is the specific heat^. The heat- 
ing rate e is given by 



(3) 



where B, and Y, are the binding energy and abundance of 
species /, is Avogadro's constant, and d/dt — S]arn{d / dx) . 
Our reaction network incorporated 430 nuclides from n to 
^'^Ge (Table[B. 

We use the reaction rates from the library R EACLIB 
dRauscher & Thielemannl2000tlSakharuk et alj2006l and ref- 
erences therein). On the timescale of the flame passage, elec- 
tron captures are unimportant, and Y^ is essentially fixed; we 
found that S]am was unchanged when weak reactions were re- 
moved from the network, so we used only strong rates for 
computational efficiency. We incorporated screening using 

^ We neglect here terms such as dE/dXj, which account for the change in 
the thermal properties as the abundance of nuclide i changes. These terms are 
much smaller than e for matter not in NSE. 



the formalism of iGraboske et alj (Il973h . Across the flame 
front, the matter does not reach nuclear statistical equilibrium 
or quasi-nuclear statistical equilibrium until most of the '^C is 
depleted. Thus, although o ur treatment of sc reening does not 
preserve detailed balance ( Calder et al. 2006) this does not af- 
fect our calculation of the flame speed. Our equation of state 
has contributions from electrons, radiation, and strongly cou- 
pled ions. We include thermal transport by both degenerate 
electrons and ph otons (for a comple te description of our ther- 
mal routines, see Bro wn et al.l l2002. and references therein). 

Equations ([T]i-(l3]l, when combined with appropriate bound- 
ary conditions, have Si^m as an eigenvalue. Ahead of the flame 
the material is at an arbitrary cold temperature Tfuei — 10^ K; 
in this region we set dT /dx to a small positive value and in- 
tegrate equations ([T]i-(l3ll. For simplicity, we split the solution 
of the thermal and network equations; that is, for each step dx 
we solve the thermal equations, to obtain T and p, integrate 
the reaction network at that T and p to compute F, and e, 
and use that to advance the solution of the thermal equations. 
For our choice of S\am, the second boundary condition is that 
F ^ Q asymptotically behind the front and that F is peaked 
where £ is maximum. We iterated until S\am had converged 
to within 0.01%. We find that S\am is insensitive to Tfuei for 
Pfuei 5 X 10^ gcm^^. At lower densities this is no longer 
true, but the relative increase in 5iam with X22 remains robust. 
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Fig. 1. — Flame speeds computed with an 130-nuclide network {dashed 
line) and a 430-nuclide networ k {dash-dotted line). We compare these 
with the results of Timmes & WooslevI 119921 dotted line), and our fit 
(eg. (4); solid line) . Our 130-nuclide network uses the same nuclides as 
ITirnines & WooslevI f 1992i) . 

Figure[T]shows a comparison between our flame sp e eds fo r 
a 1:1 C:0 mixture and those of ITimmes & WooslevI (|1992|) . 
Here we adopt the same 130-nuclides as used in that paper. 
Although in this case we are using the same nuclides and 
starting composition, the rates, equation of state, and thermal 
conductivities are not identical. Overall, our flame speeds dif- 
fer by no more than 25%; the largest discrepancy is at Pfuei = 
10^ gcm^^. Most of this discrepancy is due to different opac- 
ities used by the two codes. For Pfuei ^ 7 x 10*^ gcm^^, pho- 
tons become more efficient than electrons at transporting heat 
withi n the flame front, wit h the d ominant opacity being free- 
free. 'Tim mes & WooslevI (1199 2*) included a fit t o elec tron- 
ion scattering in the semi-degenerate regime (Ibenl l 19751) . At 
these densities where T > 2 x 10^ K and the free-free opacity 
dominates, the contribution from electron-ion scattering de- 
creases the total opacity. We compared our opacities along 
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a (p,r) trace generated for a run at Pfuei = 10 gcm^ . We 
found that non-degenerate electron-ion scattering can lower 
the opacity by « 24%, depending on how the opacity is inter- 
polated between the two limiting fits. In addition, our free- 
free opacities differ by 30% at the location along the {p,T) 
trace where \F\ is maximum. 

Finally, we also investigated the effect of reaction network 
size: increasing the network from 130 to 430 nuclides resulted 
in a 25% increase in 5iam at Pfuei = 2.0 x 10^ gcm^^; further 
increases in the size of the reaction network did not yield any 
appreciable increases in 5iani- 

3. RESULTS 

We now present the results of our flame calculations for 
different initial mixtures of '^C, '^O, and ^^Ne and different 
ambient densities. Table |2] lists S\am and the flame width de- 
fined by 5iam = {Tash - Tfusi) / max\dT /dx], with Tf^e\{ash) be- 
ing the temperature in the unburned (burned) matter. We tab- 
ulate these quantities for pg = Pfuei/10^ gcm^^ ranging from 
0.05 to 6.0, and X12 — 0.3-0.7, with the remaining composi- 
tion being '^O and ^^Ne. For each choice of pg and X12, we 
use 3 different ^^Ne abundances, X22 = 0.0, 0.02, and 0.06. 
Over most of the range in pg, X\2, and X22 in Table|2l we find 
that an increase in X22 from to 0.06 causes S\am to increase 
by approximately 30%. We confirmed several of the table en- 
tries using an i ndependent diffusion equat ion solver that uses 
adaptive grids (iTimmes & Woosleylll992h and a different re- 
action network and opacity routine. From this, we estimate 
that the uncertainty in the flame speeds listed in Table 2 are 
w 30%, with about 10% coming from numerical uncertainty 
and about 20% from different physics treatments as described 
above. We emphasize, however, that both codes find the same 
trends, e.g., an increase in S'lam with X22. 

We fit the tabulated 5iam with the approximate expression 



5iam= [23.26p9 + 37. 34p9'-^- 1.288] 
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which has fit errors, as compared against speeds calculated us- 
ing the 430-nuclide network with X22 — 0, that average 33%, 
with a maximum of 70%, for X\2 = 0.5 and pg = 0.07, and 
with a minimum of 0.1%, for X12 — 0.5 and pg — 6. For 
accurate work, interpolation from Table |2] is preferred. At 
p < 10^ gcm^^, the speedup is negligible for X12 = 0.3 but 
increases to w 50% forXi2 = 0.5. 

To understand how the addition of ^^Ne increases S'lam, we 
plot in Fig.|2]some selected abundances Y = X/A for a flame 
with an initial X12 = 0.3 and pfuei = 2.0 x 10^ gcm^^. We use 
the fraction of '^C consumed, 1 — Y12/ max(Fi2), as our co- 
ordinate and plot the region where this value is monotonic 
with distance. In a C/O deflagration, the flame speed and 
width are set by the initial burning of '^C. The buildup of 
Si-group nuclides and then establishment of nuclear statisti- 
cal equilibrium occur on longer timescales, so that the peak 
of the heat flux is reached as '^C is depleted via the reac- 
tions i2c(i2c,a)20Ne and ^^C{^^C,p)^^Na{p,af°Ne. For 
the case with X22 = 0.06 (solid lines) one sees that ^^Ne is 
depleted before the '^C is even half-consumed. The ^^Ne 
lifetime becomes less than the '^C lifetime once the a abun- 
dance is ¥4 > 10^^. At temperatures in the flame front the un- 
certainty in the ■^^Ne{a, n)^^Mg rate is estimated to be about 
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Fig. 2. — Abundances of selected nuclides during a burn at p = 2.0 x 
10' gcm^^ and with an initial '^C mass fraction of 0.3. We show runs with 
an initial ^^Ne abundance of 0.06 {solid lines) and 0.0 (dashed lines). 

10% (see iKarakas et al.ll2006l and references therein). This 
is unlike the case in AGB stars, for which the uncertainty at 
r < 3 X 10** K is approximately a factor of 10. Note that 
significant burning does not occur until the '^C lifetime be- 
comes of order the time for the flame front to move one flame 
width. This requires temperatures in excess of 2 x 10^ K 
for the densities of interest, and so ^^Ne is preferentially de- 
stroyed in a fla me via (oc,n) rather than by p-capture (cf. 
iPodsiadlowski et al. 2006). 

The neutrons made available from the destruction of ^^Ne 
capture preferentially onto ^''Ne formed during '^C burn- 
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tions build up '^O and ''^C, the latter of which then un- 



dergoes "^C(p,n)''^N(n,a)"B(p,2a)'^He. At densities of 
Pfuel ^ 5 X 10^ gcm^^ and carbon abundances X12 = 0.5, the 
flow 20Ne(n,7)2iNe(n,a)i**O(p,a)'5N(p,a)i2c also con- 
tributes. The net effect of having ^^Ne in the fuel mixture 
is that during '^C burning, the abundance of protons is de- 
pressed and the abundance of ^He elevated, as illustrated in 
Fig. ID The fact that these flows require two neutron captures 
onto the products of '^C + '^C suggests that the increase in e 
should scale roughly as Xlj. Since 5iam °= e''^^. this implies 
that the increase in flame speed will be linear in X22, which 
agrees with the numerical solution of equations ([T])-©. Be- 
cause these flows are initiated by n-capture onto ^''Ne, we 
tested our sensitivity to the reaction rate by recomputing the 
case X12 = 0.5 and Pfuei = 7.0 x 10^ gcm^^. A decrease in 
the ^''Ne + n rate by a factor of 10 produced a decrease in the 
flame speedup, from 70% to 20%. At higher densities there 
was no difference in the speedup. The only case for which 
there was no increase in 5iam was for pfuei ^10^ gcm^^ and 
X12 = 0.3 (see Table|2]i. For this case, the slower consumption 



of Ne relative to C and (n, 7) captures on Mg suppress 
the generation of a-particles early in the burn. 

4. DISCUSSION 

We have computed the laminar flame speed in an initially 
degenerate plasma consisting of '^C, '^O, and ^^Ne. We 
find that, over a wide range of initial densities and ^^C abun- 
dances, the flame speed increases roughly linearly with ^^Ne 
abundance, with the increase being « 30% forX22 — 0.06, al- 
though there are deviations from this rule at lower densities. 
These studies are relevant to the initial burning at the near- 
center of the white dwarf, and at late times where the flame 
may make a transition to distributed burning. To see how the 
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TABLE 2 
Laminar flame speed and width 
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increase in laminar flame speed changes the density where the 
burning becomes distributed, we write S'lam ~ p''(l +(^^22) 
and find from our table that at pfuei = 7.0 x 10^ gcm"^ (the 
lowest density for which our numerical scheme converged) 
and = 0.5, T? » 1.6 and t, k. 0.7/0.06. Recent numerical 
studies (IZingale et al.ll2005l) find that the Rayleigh-Taylor in- 
stability drives turbulence that obeys Kolmogorov statistics, 
so that 

"^^rn" ^i^^^^^i^^lly? wc find th3,t the fl3,nie width 
scales roughly as 5iam °= so solving for 5iam/^G = 1 im- 

plies that increasing X22 from to 0.06 would lower the transi- 
tion density by « 30%. A reduction in the density of this tran- 
sition will lead to a reduction in the mass of ^''Ni synthesized 
dHoflich et al.lll995l) . We conjecture that if a deflagration-to- 
detonation occurs, the addition of ^^Ne decreases the overall 
mass of Ni-peak elements, in addition to lowering the isotopic 
fraction of ^^Ni. 

Our results can be improved in several ways. First, the ^^Ne 
may be partially consumed as ^^C burning grad ually heats 
the core of the white dwarf (iPodsiadlowski et an[2006 ) some 



« 10^ yr prior to flame ignition. This may further reduce the 
electron fraction of the white dwarf, but will also change the 
reaction flows in the flame front. At low densities the mor- 
phology of the flame becomes more complicated, as the flows 
responsible for reaching quasi-statistical equilibrium are no 
longer fast enough to keep up with the carbon burning. In- 
deed, at p < 10^ gcm^^, the eigenfunction for the flux begins 
to show two maxima and the flame speed becomes more de- 
pendent on the ambient fuel temperature. Further studies with 
more realistic compositions and at lower ambient densities are 
ongoing and will be reported in a forthcoming publication. 
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